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In this paper we report on bubble growth rates and on the statistics of bubble topology for the 
coarsening of a dry foam contained in the narrow gap between two hemispheres. By contrast with 
coarsening in flat space, where six-sided bubbles neither grow nor shrink, we observe that six-sided 
bubbles grow with time at a rate that depends on their size. This result agrees with the modification 
to von Neumann's law predicted by J. E. Avron and D. Levine. For bubbles with a different number 
of sides, except possibly seven, there is too much noise in the growth rate data to demonstrate 
a difference with coarsening in flat space. In terms of the statistics of bubble topology, we find 
fewer 3, 4, and 5 sided bubbles, and more 6 and greater sided bubbles, in comparison with the 
stationary distribution for coarsening in fiat space. We also find good general agreement with the 
Aboav-Weaire law for the average number of sides of the neighbors of an n-sided bubble. 

PACS numbers: 82.70.Rr 



I. INTRODUCTION 

Coarsening is a process in foams by which there is diffu- 
sion of gas across films such that some bubbles grow and 
other bubbles shrink. This progresses in such a way that 
the average bubble area increases over time [1 . Coarsen- 
ing is not limited to foams, and is also relevant in other 
systems involving domain growth [3 [3] . For an ideal dry 
two dimensional foam, John von Neumann showed that 
the rate of change of area of a bubble in a two dimen- 
sional foam is 14 : 
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where n is the number of sides of a bubble, and Kq is 
a constant of proportionality. Remarkably, the shape of 
the bubble, its edge lengths, and its set of neighbors, all 
do not matter. 

In 1992 Avron and Levine [5 generalized von Neu- 
mann's law to predict the rate of area change for bubbles 
coarsening on a curved surface. The essential ingredient 
is that the sum of turning angles around each bubble is 
no longer 27r, as in flat space, but rather depends on the 
integral of Gaussian curvature, kq^ over the bubble area. 
This modifies the von Neumann law to: 



dA 



(n - 6) + - / nc^A 

7T 



(2) 



In the case of a surface of constant positive curvature, 
such as a dome of radius R, this reduces to 
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The rate of change of bubble area thus depends on the 
number of sides and the area of the bubble. 

There have been numerous theoretical and simulation 
studies of coarsening for foams in two dimensional flat 
space [6Hl4]. But to date we are aware of only one sim- 
ulation that includes the effect of substrate curvature 



[15]. There, the authors used a modified Potts model 
for two dimensional foam coarsening on spheres, toroids, 
and pseudospheres. For spheres, they focused on how 
the area distribution and average area change over time, 
and find that at late times the dynamics are dominated 
by the appearance of 'singular bubbles' much larger than 
the average that quickly grow to cover the sphere. There 
is minimal discussion of the coarsening of individual bub- 
bles, and no discussion of side number or other distribu- 
tions of the system. 

While the coarsening of foams in two dimensional flat 
space has been well-measured [12 , 16 - 25 , we are unaware 
of any experiments to test the modified law of Avron and 
Levine for foam in two dimensional curved space. How- 
ever metallurgical grain growth on curved substrates has 
been reported. In Ref. [26 , the results are said to be pre- 
liminary and no growth rate data are shown. In Ref. [27] , 
the deviation from the coarsening rate for flat space is 
masked by noise, but statistical analysis is reported to 
demonstrate consistency with Eq. In this paper we 
use a hemispheric cell to create a curved two dimensional 
foam. We use image analysis to track individual bubbles 
and measure bubble dynamics such as coarsening rate. 
We also measure bubble statistics, such as the distribu- 
tion of number of sides and compare this to results from 
a fiat cell. Our image quality and analysis methods are 
sufficient to demonstrate directly, for six-sided bubbles, 
that the growth rates are different from flat space and 
are consistent with Avron and Levine [5]. 



II. MATERIALS AND METHODS 

To measure coarsening rates of two dimensional foams 
on a curved surface, we constructed a cell from two hemi- 
spherical polycarbonate domes. The smaller dome has an 
outer diameter of 12.5 cm, and the larger dome has in- 
ner diameter 13.3 cm, creating a 4 mm gap. The smaller 
dome was glued to a flat acrylic plate. The larger dome 
was placed over the smaller dome and separated from 
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FIG. 1: Sample photograph of a two-dimensional foam, coars- 
ening between nested polycarbonate hemispheres with a 4 mm 
gap. The scale bar is 1 cm. 

the plate by an 0-ring of diameter 0.25 inches. We were 
careful to ensure that the two domes were aligned con- 
centrically. The upper dome was then screwed to the 
plate to create a sealed chamber of constant curvature. 

The solution we used to create our foam was a liquid 
consisting of 75% deionized water, 20% glycerin, and 5% 
Dawn Ultra Concentrated dishwashing liquid. This cre- 
ated a foam that was stable and generally lasted many 
days. The foam was prepared by putting 35 mL of solu- 
tion into the chamber (this fills the dome to about 2 cm 
above the 0-ring) and shaking it until a uniform opaque 
foam was created, with an average bubble size much less 
than the separation of the domes. The chamber was then 
left to coarsen until a single layer of bubbles remained be- 
tween the two domes. This took about 24 hours. Two 
dimensional coarsening could then be observed for the 
next 2-4 days. Film ruptures were sometimes observed 
at the end of this period. Bubble statistics were not taken 
after any ruptures, though single bubble dynamics were 
still considered. 

To photograph the foam, the chamber was placed 
65 cm above a Vista Point A lightbox, which provided a 
spatially and temporally uniform light source. A Nikon 
D80 camera with a Nikkor AF-S 300 mm 1:2.8 D lens 
was mounted 2.5 m above the chamber. The camera was 
controlled by a computer to take pictures every two min- 
utes. The apparatus was left to collect pictures, for a 
period ranging from a few days up to a week. This pro- 
cess was repeated three times to build up statistics. A 
sample photograph can be seen in Fig. [l] Note in this 
photograph that, especially towards the edge, it is possi- 
ble to see both the Plateau borders on the top dome, as 
well as the Plateau borders on the bottom dome. This 
makes it difficult to identify the correct boundaries of the 
cells. We address this issue as part of our image analysis. 

After we have taken a series of images of the dome, we 
perform analysis to get out the true areas of the bubbles 
on the dome. Note that the images constitute an or- 
thographic projection of a sphere (or hemisphere) onto a 
plane, where the point of projection is infinity. This pro- 



FIG. 2: Transformation of the photograph from Fig. ^ by the 
Lambert equal area projection, to get the true areas of the 
cells. There are errors in this result, but for a region in the 
center, typically about 4.5 x 4.5 cm^, it works well and the 
areas of these cells can be tracked. This image is skeletonized 
before any actual areas are measured. The scale bar is 1 cm. 



jection converts the positions according to the following 
equations [28^: 

X = R cos (p sin A , x 

y = R sin cp ^ ^ 

where R is the radius of the sphere, A is the longitude, (f 
is the latitude, and the center of the domes is defined as 
X = 0,(p = 0. The first problem is that in a given image, 
both the Plateau borders on the top and bottom domes 
are visible. In order to isolate a single set of Plateau 
borders so that the cells' edges are defined correctly, we 
recognize that the image is projected in two ways. The 
Plateau borders on the top dome are an orthographic 
projection of the foam using the radius of the top dome, 
and the Plateau borders on the bottom dome are an or- 
thographic projection of the foam using the radius of the 
bottom dome; both are combined into the same image. 
To undo this transformation, we take the inverse of the 
transformation twice, once using the radius for the top 
dome and once using the radius for the bottom dome. 
The resulting two images are thresholded and dilated. 
The images are then multiplied. This kills the Plateau 
borders that do not correspond to the transformation. 
That is, the Plateau borders from the bottom dome that 
were transformed using the radius of the top dome are 
killed and vice versa. The result is a binary image with 
the correct latitudes and longitudes of the Plateau bor- 
ders on the dome. 

After we have accounted for the fact that the Plateau 
borders on both the inner and outer domes are visible, 
we can then consider the areas of the individual bubbles. 
The binary image with the correct latitudes and longi- 
tudes of the Plateau borders has errors, but does well 
for a region of interest in the center. This resulting im- 
age, however, does not preserve the areas of the cells. 
A simple projection that will preserve areas is the Lam- 
bert cylindrical equal area projection. This projection is 
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FIG. 3: (Color online) Area versus time for three different six- 
sided bubbles, with initial area, Aq, as labeled. Lines are a 
linear fit to the data, from which Aq was found and subtracted 
from the data. The vertical error bars represent perimeter 
times pixel size divided by the square root of number of points 
in the perimeter; this statistical uncertainty matches the scat- 
ter in the data points. 



defined by [28]: 



x = RX 
y = R sin (p 



(5) 



where as before A is the longitude and cp is the latitude. 
The R used here is the average of the two domes used. 
This produces an image with cells that have distorted 
shapes, but have the same areas as the actual cells on 
the curved surface. This allows us to track individual 
bubble areas over time. The result of this image analysis 
can be seen in Fig. |2] 



III. BUBBLE DYNAMICS 

Using the method of finding the correct areas of in- 
dividual bubbles described above, it is possible to track 
the areas of individual bubbles over time. The method 
of identifying the correct Plateau borders and finding the 
correct areas sometimes has errors of failing to identify 
a Plateau border or adding an extra one, especially far- 
ther from the center, where distortion is greater. It was 
possible to find bubbles, especially near the center, that 
would be correctly identified for a significant length of 
time. The viewing region where bubbles can be tracked 
is typically about 4.5 x 4.5 cm^. Not all bubbles in this 
region will necessarily be able to be tracked, but it is very 
rare for a bubble outside this region to be trackable. The 
projections were compared to the original images, to en- 
sure that the areas of the tracked bubbles were correct. 
Images were taken at 200 minute intervals and analyzed 
to get the correct bubble areas. Correctly identified bub- 
bles were labeled and tracked over this period. From this 
it was possible to get area versus time graphs for many 
bubbles. 

We can now consider the form that the area versus time 
traces for the bubbles should take. Avron and Levine's 
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FIG. 4: (Color online) Growth rate of bubble area versus the 
expected factor from Avron and Levine, Eq. (|3|, where A is 
bubble area and R is the radius of the dome. The error bars 
represent the uncertainty in slopes of linear fits to bubble area 
versus time. The solid line is a proportionality fit to all data, 
with slope Ko = 0.00018 ±0.00008 cm^min. The inset shows 
a blow-up of the data for the n = 6 sided bubbles. There, the 
growth rates for the three representative bubbles of Fig. [S] are 
circled in color. 



prediction for the coarsening of foam on a spherical sub- 
strate of radius R^ Eq. (3|, is a linear differential equation 
that can be solved analytically for area versus time: 



Ao + Kot[n-Q 
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Here Ao is the area of a bubble at time zero, and Eq. ^ 
is the Taylor expansion in {Kot/ R^). Note that the stan- 
dard von Neumann result, A = Aq -\- Ko{n — 6)t, is ob- 
tained in the limit R ^ oo. 

Examples data for the area versus time of three six- 
sided bubbles with different initial areas are plotted in 
Fig. |3] There the initial area of each bubble was sub- 
tracted off so that the traces are easily comparable. The 
lines are a linear fit to the data, giving a constant growth 
rate dA/dt that is positive. It is possible to fit the data 
to the full exponential form of Eq. ([6|, but the addi- 
tional terms in the expansion from Eq. ([7| are much, 
much smaller than the linear term, so it is sufficient to 
fit to an ordinary line. The key features of Fig. [3] is that 
the six-sided bubbles grow, and that the larger ones grow 
faster. This agrees with Avron and Levine, and contrasts 
strongly with the case of a flat sided cell, where six-sided 
bubbles neither grow nor shrink according to the usual 
von Neumann equation. In particular, in our recent ex- 
periments on the coarsening of bubbles in a flat cell [25] , 
where the liquid fraction was varied, the six-sided bub- 
bles all had growth rates scattered around zero to within 
statistical uncertainty. 

We now measure the growth rate for all bubbles, as 
illustrated in Fig. [3) and we compare to the expected 
relationship from Avron and Levine's modification to 
von Neumann's law. This is plotted in Fig. |4j where 
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FIG. 5: (Color online) Growth rate of bubbles versus area. 
Symbol types distinguish bubbles with different number n 
of sides, as labeled. Horizontal lines are the expected rela- 
tionship from the uncorrected von Neumann law, shown in 
Eq. ([T]), using Ko from Fig. [4] Solid lines are the expected re- 
lationship using the corrected von Neumann law from Avron 
and Levine, shown in Eq. (|3|, using Ko from Fig. [4] The data 
and vertical error bars are the same as shown in Fig. ^ The 
horizontal error bars represent the range of area over which 
each bubble grew with the specified rate and number of sides. 



the y-axis is the coarsening rate, and the x-axis is the ex- 
pected proportionality for a dome of constant curvature 
given by Eq. (|3|. Each point represents a single bubble. 
The line is a proportionality, with slope Kq, which is the 
only fitting parameter. The inset is a blow-up showing all 
the six-sided bubbles, and highlighting the three bubbles 
featured in Fig. [3j For six-sided bubbles the growth rates 
are all positive, except for one or two outliers. There is 
notable scatter, but the evident trend is that dA/dt in- 
creases with bubble size. 



Another way to compare growth rate data to Avron 
and Levine is to plot the coarsening rate against area, as 
shown in Fig. [5] Here each point represents a single bub- 
ble, color coded by the number of sides. The horizontal 
lines are the expected relationship from the unmodified 
von Neumann's law, as seen in Eq. 0, using; the same 
constant of proportionality Ko as found in Fig. [4j The 
solid lines are the expected relationship from the modi- 
fied von Neumann's law, as seen in Eq. (|3|, agam usmg 
the same value of Kq. We see that the data are generally 
consistent with the predicted modification. This is most 
evident in the six-sided bubbles, which are virtually all 
growing. The rate of area change also appears to increase 
with area for n = 7, but is masked by noise for other n. 



IV. BUBBLE DISTRIBUTIONS 

With our system it was also possible to measure dis- 
tributions such as p(n), the probability that a bubble 
has n sides, and m{n)^ the average number of sides of 
the neighbors of an n-sided bubble. Unlike the case of 
the flat cell, we do not expect to reach a scaling state 
where these statistical quantities remain constant over 
time. Because the growth rate of a bubble grows with its 
area, we expect at long times to have large bubbles grow 
rapidly to dominate the system, and this will cause bub- 
ble statistics and distributions to change with time. Our 
system is at much earlier times, where the modification 
to a bubble's growth rate due to its area is small. This 
modification still should have some impact and we do in- 
deed find that our statistics deviate from the ordinary 
scaling state observed in the fiat cell. 

To measure our statistics, a 8 cm x 8 cm region of 
interest was defined in the center of the dome, and the 
number of sides of all bubbles completely or partially 
within this region was recorded by hand. This was done 
for photographs at 400 minute intervals from the earliest 
photograph of a two dimensional foam to the first occur- 
rence of a rupture, for a total period of typically 48 hours. 
This was repeated for three runs. This data was used to 
produce a distribution of the number of sides, which can 
be seen in Fig.[6j Also shown for comparison is the distri- 
bution found for a flat two dimensional foam in Ref. [25] . 
We see that as compared to the flat cell, the dome has a 
surplus of six-sided bubbles, and a deficit of 3, 4, and 5 
sided bubbles. 

It is also possible to describe these distributions by 
their average, (n) = 6.06 ±0.1, and by their variance, 
/i2 = — 6)^. We measure the variance to be 

/i2 = 1.30 ± 0.05. This value is lower than was measured 
for the flat cell, fi2 = 1.56 ±0.02 [25 , indicating that the 
width of the distribution is narrower. 

From this same data we can measure m(n), the average 
number of sides of the neighbors of an n sided bubble. 
We expect m{n) to be related to n by the Aboav-Weaire 
law, which predicts m{n) = 6 — a ± [(6a — /i2)/^]- In 
this equation a is a fitting parameter which is usually 
found to be around 1. Our measurements for m(n) can 
be seen in Fig. [t] along with the measurements of m(n) 
for a flat cell. Fits to the Aboav-Weaire law are also 
shown, using the relevant value for /i2 in each case. We 
find a = 1.1 ±0.08 for the flat cell [25^ and a = 0.96±0.09 
for the dome. We see that in both cases there appears to 
be more curvature in the data than predicted. The data 
for the flat cell also seems to fit the form better than for 
the dome. 



V. CONCLUSION 

In this experiment we were able to measure both 
bubble statistics and bubble dynamics of a foam on a 
curved two-dimensional surface of radius R. For bubble 
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statistics we find that bubbles with few sides are under- 
represented as compared to a two dimensional foam in 
flat space. We also find that the Aboav-Weaire law gen- 
erally holds, though not quite as well for the dome as 
for the flat cell. For bubble growth rates, in general, 
it is difficult to observe the effect of the term added to 
von Neumann's law by Avron and Levine to account for 
substrate curvature. Our maximum bubble size is around 
^max = 3.5 cm^, as compared to = 41.6 cm^; there- 
fore, for all our bubbles 3A/{7rR'^) <C \n — 6\ holds, except 
for n = 6. This is why all the data in Fig.[4]lie at x values 
near (n — 6) . Even if we managed to get a single bubble 
of 20 cm^ to completely fill our viewing area, the rate of 
area change would be dA/dt = Ko[{n — 6)-\-0A6], so that 
a seven sided bubble would have less than a 50% increase 
in growth rate. For these reasons, the clearest signal we 
see of curvature effects is that six-sided bubbles system- 
atically grow and do so faster for larger bubbles. The 
coarsening data as a whole is consistent with Avron and 
Levine's modification to von Neumann's law to account 
for coarsening on a curved surface. 



FIG. 6: (Color online) Side number distribution for a flat cell 
and for the dome. For the flat cell, the average number of 
sides is (n) = 5.92 ± 0.01 and the variance /i2 = 1.56 =b 0.02. 
For the dome (n) = 6.06 ±0.1 and /i2 = 1.30 ±0.05. Data for 
flat cell is taken from Ref. J25^. The vertical error bars are 
from a fractional area of one over the square-root of number 
of bubbles with specified side number. 
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FIG. 7: (Color online) Aboav-Weaire relationship for a flat 
cell and for a dome. The y-axis is m(n), the average number 
of sides of an n sided bubble. Lines are fits to the Aboav- 
Weaire law, m(n) = 6 — a ± [(6a — iJ.2)/n], where a is the only 
fitting parameter. We find a = 1.1 ± 0.08 for the flat cell and 
a = 0.96 ±0.09 for the dome [25 . The variances, /i2, for these 
systems can be found in the caption to Fig. [6] 
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